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Abstract: In Raman spectroscopy of graphite and graphene, the D band at ~ 1355cm -1 is 
used as the indication of the dirtiness of a sample. However, our analysis suggests that the 
physics behind the D band is closely related to a very clear idea for describing a molecule, 
namely bonding and antibonding orbitals in graphene. In this paper, we review our recent 
work on the mechanism for activating the D band at a graphene edge. 
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;> ! 1. Introduction 

X- 

Bonding and antibonding orbitals are basic ideas for describing molecules. Bonding orbitals 
contribute to the formation of a molecule, whereas antibonding orbitals weaken the bonding and 
destabilize a molecule. Normally, bonding orbitals are more stable than antibonding orbitals in terms of 
energy and thus a molecule is stable unless sufficient electrons occupy the antibonding orbitals. 

Graphene [1,2] is unique with respect to its molecular orbitals. The bonding and antibonding orbitals 
in graphene are degenerate, and various types of linear combination of these orbitals form the Fermi 
surface expressed by the isoenergy sections of Dirac cones. This degeneracy plays an essential role 
in various phenomena. For example, graphene is stable with respect to a large shift of the Fermi 
energy position [3]. Another notable example is that graphene exhibits high mobility. Elastic backward 
scattering between the bonding and antibonding orbitals induced by long-range impurity potential is 
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suppressed because they are orthogonal [4]. In this paper, we show that the bonding and antibonding 
orbitals in graphene are key factors in the activation mechanism of the D band observed at a graphene 
edge. 

Since the discovery of the D band, much interest has focused on its origin. Tuinstra and Koenig 
attributed the D band to an A lg zone-boundary mode at the edge of a sample (see Fig. 1) on the grounds 
that the Raman intensity is proportional to the edge percentage and that the edge causes a relaxation of 
the momentum conservation needed for activating a zone-boundary phonon [5]. Katagiri et al. confirmed 
that the D band originates from an edge (or discontinuity in the carbon network) by observing the light 
polarization dependence of the D band intensity at graphite edge planes [6] . The atomic arrangement 
of an edge has two principal axes; armchair and zigzag edges. Cancado et al. showed that the armchair 
(zigzag) edge is relevant (irrelevant) to the relaxation of the momentum conservation for a zone-boundary 
phonon [7] . In addition, they found that the D band Raman intensity depends on the polarization of laser 
light, that is, the intensity is maximum (minimum) when the polarization is parallel (perpendicular) to 
the armchair edge. The light polarization dependence of the D band is also observed ubiquitously at the 
armchair edges of a single layer of graphene, which suggests that out-of-plane coupling in graphite is 
not essential to the origin of the D band [8-10]. 

A model of the D band must at least explain the observed properties: (i) the D band intensity increases 
only at an armchair edge, and (ii) the D band intensity is dependent on the laser light polarization. 

The current D band model is a double resonance model [11]. In this model, a photo-excited electron 
passes through two resonance states, which enhances the Raman intensity of a phonon with nonzero 
wave vector q ^ 0. This model is not concerned with the details of electron-phonon and electron-light 
matrix elements, and it does not provide clear explanations of the properties of the D band. Also, the 
intensity calculated with this model is dependent on the lifetime of the resonance states. Usually, the 
lifetime is determined in such a manner that a calculated result reproduces experimental data. In this 
sense, the double resonance model is phenomenological. Because the lifetime is shorter in a defective 
graphene sample, double resonance does not necessarily mean an enhancement of the D band Raman 
intensity. On the other hand, the model can account for the so-called dispersive behavior of the D 
band [11]. However, as we will show in this paper, dispersive behavior is characteristic of A lg modes, 
rather than an inherent property of the D band. In fact, dispersive behavior is observed also for the 2D 
band, and the excitation does not need an edge, which is in contrast to the D band. 

In this paper we show that the observed properties of the D band are naturally explained in terms of 
simple ideas based on molecular orbitals and momentum conservation. In our formulation, the D band 
is excited from a photo-excited electron through a single resonance process in the same way as the G 
band. It is concluded that, without invoking an artificial assumption, the D band is closely related to the 
orbital dependence of the electron-phonon matrix element, the special nature of the armchair edge, and 
optical anisotropy [12]. 

This paper is organized as follows. In Sec. 2 we observe that a graphene molecular orbital and wave 
vector are closely correlated. This correlation is both an important factor in terms of understanding the D 
band and an essential feature of graphene. In Sec. 3 the properties of the D band are deduced from three 
factors. The importance of the electron-phonon matrix element and the role of the armchair edge in the 
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excitation mechanism of the D band are explained in detail. In Sec. 4 we show some predictions obtained 
with our model. Future prospects and our conclusion are given in Sec. 5 and Sec. 6, respectively. 

2. Bonding and antibonding orbitals 

Graphene's hexagonal unit cell has two carbon atoms, denoted by A and B in Fig. 1, and the electron's 
wave function ip is written as a linear combination of 2p z atomic orbitals of the A and B atoms, xa and 
Xb- When we apply Bloch's theorem to graphene, we obtain tp and the band structure as a function of 
the wave vector k [13,14]. In the Brillouin zone (BZ) of graphene, there are two points, namely the K 
and K 7 points, where the conduction and valence bands touch each other. The orbitals of the states near 
the K point take the form of 



where the phase e [0, 2ti] is the polar angle between vector k measured from the K point and the 
/c x -axis (see Fig. 1), and s = ±1 is the band index (s = +1 is the 7r*-band and s = —1 7r-band). The 
orbitals of the states near the K' point are written as 



where Q' is the polar angle defined with respect to the K 7 point as shown in Fig. 1. 

Figure 1. Graphene unit cell and BZ. The K (K 7 ) point is located at k F = (47r/3a, 0) (— k F ), 
where a is the lattice constant. The positions of the bonding and antibonding orbitals are 
marked by red and blue points on the iso-energy sections of the Dirac cones (circles). The 
Raman D band is composed of zone-boundary A Xg modes that consist only of C-C bond 
stretching motions. The armchair edge identifies a state at (k x , k y ) with a state at (— k x , k y ), 
and causes the BZ holding. 



^(k) = i=( e - e « XA + SXB ), 
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In Eq. (1), the bonding and antibonding orbitals (xa + Xb and — Xa + Xb) are located at = and n, 
respectively, on the iso-energy section of the 7r*-band (s = +1). The orbital with a general is a linear 
combination of the bonding and antibonding orbitals. In Eq. (2), the bonding and antibonding orbitals 
are located at 0' = n and 0, respectively, on the iso-energy section of the 7r*-band. As we can see in 
Fig. 1, the bonding (antibonding) orbitals are located symmetrically on the /c^-axis with respect to the Y 
point, which is one of the most important characteristics of the BZ of graphene. 

The orbitals Eqs. (1) and (2) can be derived from an effective model for graphene, that is, a massless 
Dirac equation, in the following manner. The energy eigenequation is written as 



Ei; = Hi; = v F \ * 1 , + , (3) 




where v F is the Fermi velocity, p = (p x ,p y ) is a momentum operator, in the off-diagonal terms 
represents a 2 x 2 null matrix, and cr = (a x , a y ) and cr' = (—a x , a y ) are 2 x 2 Pauli spin matrices: 





o)- <4) 

The massless Dirac equation in Eq. (3) is decomposed into two Weyl's equations for the K and K 7 
valleys: Eip K (r) = v F er ■ ptp K (r) and E^ K i(r) = v f <t ' • p^K'( r )' respectively. To reproduce Eq. (1), 
we assume the plane wave solution ^(r) = e lk ' T i/j K (k)/y/V, and obtain the energy eigenequation 
hvpcr ■ k0 K (k) = Et/j K (k). Then, the corresponding energy eigenvalue and eigenstate are easy to find 
using k x =f ik y = |k|e Ti0 ^ as E = shv F \k.\ and 

i / p -i0(k)\ 

This is identical to Eq. (1) by setting 

xa=Q, »=(?)■ 

Similarly, we can check that Eq. (2) is the solution of hv F er' ■ k0 K '(k) = Eip K >(k), which is given by 

1 /_p*e'(k)\ 

*«^( s J- (7) 

The Dirac equation is very helpful as regards understanding the mechanism of a result in terms of 
symmetry and momentum conservation. In particular, the fact that the Dirac equation is composed of a 
multiplication of the Pauli matrices (for the A and B atoms) and momentum operators makes it easy to 
recognize that the orbital (the pattern of the linear combination of xa and xb) is dependent on the wave 
vector of the particle. In the following, we refer to the Dirac equation for graphene in order to capture 
the essential features of a result. The graphene Dirac equation differs from the original Dirac equation in 
the following respect: the wave function ip in the graphene Dirac equation has 4 components consisting 
of two orbitals (xa and xb) and two valleys (K and K'), while those in the original Dirac equation are 
two spin states (up and down spins) and two chiralities (left- and right-handed) [15]. Thus, the orbital 
degree of freedom is commonly referred to as the pseudo-spin. 
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3. Light polarization dependence of D band intensity 

In this section we show that the light polarization dependence of the D band originates from three 
factors. The first factor concerns the nature of the electron-phonon interaction for the A lg mode, which 
will be explained in Sec. 3.1. The second factor is a modification of the BZ by the armchair edge, which 
will be explored in Sec. 3.2. In Sec. 3.3, we describe the third factor, which concerns the interaction 
between electrons and a polarized laser light. In Sec. 3.4, we construct the D band polarization formula, 
by combining the three factors. 

3.1. Dominance of intervalley backward scattering 

Suppose that an electron has been excited into the 7r*-band by a laser light [Fig. 2(a)]. When 
a photo-excited electron emits an A ig mode, there is a strong probability that the electron will 
undergo (intervalley) backward scattering, as shown by the real space diagram in Fig. 2(b). In 
the A;-space, the change in the exact (approximate) intervalley backward scattering is denoted by 
the black (orange and green) solid arrow. Although the forward scattering denoted by the dashed 
arrow may be allowed by momentum conservation, it never takes place because orbitals suppress the 
corresponding electron-phonon matrix element. This dominance of intervalley backward scattering 
originates from the characteristic feature of an A\ g mode, namely that the vibration consists only of 
bond shrinking/stretching motions, as shown by the displacement vectors in Fig. 1. Mathematically, this 
characteristic of the A lg mode is described by the fact that the electron-phonon interaction, H ep (D), 
satisfies 

XaH cp (D)xb = XbH cp (D)xa = 9e P , (8) 
X {H cp (D)xa = xiH cp (D) X B = 0, (9) 

where g ep is a coupling constant for bond stretching. With these H ep (D) conditions, we obtain the 
electron-phonon matrix element squared \M\ 2 = \^,(Q'y H ep (D)ijj^(Q)\ 2 using Eqs. (1) and (2) as 

\m\ 2 = ^{i- ss 'cos(e'-e)}. (io) 

Equation (10) shows that, for intraband scattering (ss' = 1), the scattering probability of the exact 
backward scattering (0' = + n) is maximum, while that of the exact forward scattering (©' = 0) 
vanishes. 

When the energy of a photo-excited electron is much larger than the phonon energy (~ 0.15 eV), we 
can assume that no significant energy shift occurs as a result of the inelastic scattering. For the exact 
backward scattering, the A lg wave vector relates to the photo-excited electron wave vector as q = —2k, 
where q (k) is the wave vector of the A lg mode (photo-excited electron) measured from the K point. 
This relationship between k and q shows that the orbitals effectively relate the electron wave vector k to 
the Aig wave vector q. 1 

'Thus, the orbital dependence of the electron-phonon matrix element justifies the basic idea of the quasi-selection rule 
for the D band proposed by several authors [16,17]. As examined in Sec. 4.1, the relationship between k and q ( q = —2k) 
is important when we discuss the dispersive behavior of the D band. 2ftt>F|k| is approximately equal to the laser excitation 
energy El and fi^F|q| changes linearly with changing El- 
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Figure 2. Mechanism of D band intensity enhancement at armchair edge, (a) A 
photo-excited electron near the K point, (b) The electron-phonon interaction of an 
zone-boundary A\ 3 mode results in the dominance of intervalley backward scattering, (c) 
The BZ holding leads to the appearance of two special electronic states (with the bonding 
[red] and antibonding [blue] orbitals): an A\ g mode can be excited from them through a 
first-order Raman process. In real space, this is represented by the fact that the color of the 
points does not change while emitting an A lg mode. 
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It is straightforward to reproduce Eq. (10) in the framework of the graphene Dirac equation. The 
electron-phonon interaction for an A lg mode with wave vector q is written as [18] 

(0 e~ iqr (T A 

and the matrix element M is given by 

M= f rf 2 r^',(r)ti7 cp ( J D q )^(r) 
Jv 




The last line is written as a multiplication of two parts: the first part represents momentum conservation 
and the wave vector of the scattered electron k' is given by k' = k + q. The second part gives Eq. (10). 
In addition to momentum conservation, we can use energy conservation to obtain i7 F jk'| = i7 F |k| — w q , 
where w q is the frequency of A lg . Thus, when ftu F |k| 3> ftw v we have |k'| ~ |k|. Since the orbital part 
results in the dominance of intervalley backward scattering, we obtain k' ~ — k. As a result, q ~ —2k 
is satisfied. 

3.2. Brillouin zone holding 

The dominance of intervalley backward scattering causes an enhancement of the D band Raman 
intensity if Brillouin zone holding (BZH) by the armchair edge is taken into account [12]. Here, BZH 
means that, in the BZ of graphene shown in Fig. 1, a state with (k x , k y ) is identical to a state with 
(— k x , k y ) and that the correct BZ is given by the positive k x region of the original BZ of graphene [19]. 
In Fig. 2(c), as a consequence of BZH, the final state in the intervalley backward scattering event is 
identified with the state near the K point. Namely, the state with (— \k x \, —k y ) near the K 7 point is 
identified with the state with {\k x \, —k y ) near the K point. So, in the held BZ, the change in a photo- 
excited electron is (\k x \, k y ) —¥ {\k x \, —k y ), as shown by the solid arrow in Fig. 2(c). Generally, the 
probability of a process in the held BZ is given by replacing 0' with it — 0' in Eq. (10) as 

Q 2 

|M B zh| 2 = ~f {l + ss'cos(e' + e)}, (13) 

and this transition probability is indeed maximum when ©' = —0. When 0' = and s' = s in Eq. (13), 
the final state coincides with the initial state, which is the condition of a first-order Raman process. In 
this case, the probability is given by 

|M B zh| 2 = ^ p cos 2 0, (14) 

which takes its maximum value for = and it. This shows that the states near the /c^-axis (or the 
states with bonding and antibonding orbitals) can contribute to the D band intensity through a first-order 
Raman process. 

BZH originates from the fact that a special standing wave is formed by an armchair edge. The standing 
wave is constructed by the antisymmetric combination of an incident plane wave with k = (k x , k y ) and 
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a scattered plane wave with k' = (- 



■ k k ~\ 

n X) n y ) 



as e 



oc e lkyV sin(k x x). Asymmetric 



combination does not satisfy the boundary condition for the armchair edge. Note that sm(k x x) does 
not change when k x is replaced with — k x , except for the unimportant change in the overall sign. More 
importantly, the orbital part Eq. (1) does not change when there is the reflection at the armchair edge, as 
we can confirm by replacing 6' with ix — 6 in Eq. (2) [20]. The same orbitals are superposed to form a 
standing wave at the armchair edge. Thus, the total wave function becomes 



V£(k) = N < e 



ikr 



-i9(k) 



ik'-r 



N > e ik y y sjn^z) 



-*9(k) 



(15) 



where N (N') is a normalization constant [19]. The standing wave does not change with the replacement 
k x — k x , and therefore the correct BZ of the standing wave is given solely by the positive k x region to 
avoid double counting. It is noteworthy that BZH is specific to the armchair edge and is not applicable to 
a zigzag edge. The absence of BZH at a zigzag edge is due to the orbital part changes with the reflection 
of an electron at the zigzag edge and the different orbitals are superposed to form a standing wave. The 
standing wave for a zigzag edge is written as 



^(k) = Ne ik * x < e ik * y 



-i9(k) 



N'2ie 



/ I\ T J. 



sm(k y y - ©)) 
ssin(k y y) 



(16) 



which is not invariant with the replacement k y — > —k y unless = or tt. Thus, it needs a second-order 
process to activate a phonon mode with nonzero q (see Fig. 3) and a first-order Raman band (except the 
G band) cannot appear at the zigzag edge. Although the intensity is not comparable to that of G or D 
bands, we can expect a second-order band to appear at the zigzag edge. The D' band [21,22] (not the D 
band) may be such a second-order Raman band that can be described by the double resonance model. 

Figure 3. A process that causes the D' band at the zigzag edge. Phonon excitation (q ^ 0) 
accompanied by a change in the electronic states, which shows that the phonon does not 
appear through a first-order process. Note that, in the right diagram, the color of the points 
changes after the phonon is emitted. 




Zigzag Edge 
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The standing wave at the armchair edge can be expressed in the framework of the Dirac equation as 

/ e +ik x x g-i0(k)\ 



< k (r) = Ne*** 



^-\-ik x x g 
e -ik x x e ~i&(k) 




, p +ik x x\ i /„-ie(k) 



where k denotes the wave vector measured from the Dirac point in the held BZ and ® represents the 
direct product of the valley and the orbital. It can be confirmed that Eq. (17) reproduces Eq. (13), by 
calculating the expectation value of i7 e p(-Dq) of Eq. (11) with respect to Eq. (17), 

M BZH = / d 2 r< k ,(r)t4 p (7J q )^ k (r) 
Jv 

= 5(q x -k x - k' x )\ {5(k y + q y - k' y ) + 5(k y - qy -k' y )} 9 -f (e l& s + s'e^) , (18) 

where we used N 2 f dye^^'^ f dxe^--^'^ = (l/2)6(k y ± q y - k' y )S(q x - k x - k' x ). For 
k' = k (i.e., for a first order Raman process), Eq. (14) is reproduced, and momentum conservation gives 
q=(2fc x ,0). 

3.3. Optical anisotropy 

The third factor is the polarization dependence of optical transitions. To supply photo-excited 
electrons to the states with bonding and antibonding orbitals on the /c^-axis, the polarization of incident 
laser light must be set parallel to the armchair edge (y-axis), as shown in Fig. 4(b). The photo-excited 
electrons on the A^-axis emit A lg modes without changing their positions [see Fig. 4(a)]. Meanwhile, 
when the polarization of the incident laser light is set perpendicular to the armchair edge, the rr-polarized 
light supplies the states on the A; y -axis with photo-excited electrons [see Fig. 4(c)]. The electrons near 
the fc y -axis change their positions in the held BZ when they emit A lg modes [see Fig. 2(b)], and these 
electrons on the k y -axis do not contribute to the D band intensity. As a result, the D band can be 
strongly dependent on the laser light polarization: the D band intensity is enhanced (suppressed) when 
the polarization of the incident laser light is parallel (perpendicular) to the armchair edge. 

The optical anisotropy is due to the dependence of the optical matrix element [23]. Namely, for y 
(x) polarized light, E y (E x ), the optical matrix element squared is proportional to cos 2 6 (sin 2 6) as 

|M opt | 2 =( ^ C0S2@ f ° r4 ' (19) 
# 2 7 sin 2 e for 7^, 

where g cl denotes an electron-light coupling constant. The electrons near the k x (fc 2/ )-axis are dominantly 

— * — * 

photo-excited by E y (E x ) [see Fig. 4(b) and (c)]. For the general polarization direction of incident laser 
light E, M opt is proportional to the direct product of E and k (= k/|k|) as M opt ocExk. 

Although Eq. (19) was derived without taking account of the edge, it turns out that a similar optical 
matrix element is obtained for the standing waves [24]. Here, let us use the Dirac equation of Eq. (3) to 
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Figure 4. Light polarization dependence of the D band, (a) An A\ g phonon is excited 
through a first-order process from the electrons with bonding and antibonding orbitals on 
the fc x -axis. (b) To supply a photo-excited electron to the states on the k x axis, the laser 
light polarization must be parallel to the y-axis or the armchair edge, (c) E x can produce a 
photo-excited electron at the A^-axis. Such an electron changes its position when it emits an 
Ai g mode, and thus it cannot contribute to the D band. 




obtain the matrix element that includes the effect of the armchair edge. The electron-light interaction is 
given by replacing p with p — eA in Eq. (3) as 

y a' ■ (p - eA)) 

where A is the vector potential of light. Since the Maxwell equation gives E = —dA/dt, the vector 
directions of E and A are the same. The optical matrix element is defined using Eq. (17) as 

M opt = ~ev F J dh^(ry ^ q A J \ V^(r). (21) 

For the y-polarized light A = (0, A y ) , M opt is nonzero only for a direct transition (k' = k) and the orbital 
gives cos6. For the x-polarized light A = (A x ,0), M opt includes the integral (1/L) dxsm(k x — k' x )x 
and the orbital gives (— e* ' + e~ l6> )/2. Since the integral vanishes when k x = k' x , a direct transition 
does not take place. The possible transitions are indirect transitions k x ^ k' x . For k x — k' x = rnr/L (n 
is an odd number), we obtain 2 /(nit) by performing the integral. The momentum change in an indirect 
transition is inversely proportional to the distance (L) from the armchair edge where the standing wave 
is a good approximation. Since the change in the wave vector is negligible when L is large, we may 
assume 6' ~ 0. Then the orbital leads to — ism 6, which reproduces Eq. (19). This feature of the 
indirect transition for the x-polarized light has been examined with a more mathematically rigorous 
method using a lattice tight-binding model [24]. 



3.4. D band polarization formula 



We combine Eqs. (14) and (19) to derive the polarization formula of the D band. The probability 
of a first-order Raman process that an electron with 6 in the 7r-band is excited into the 7r*-band by E { 
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(i = {x, y}), and then the photo-excited electron emits the A lg mode, and finally the electron with 
emits a light with Ej is given by 

{cos 6 6 (ji) = {yy) 

cos 4 9 sin 2 9 (xy) or (yx) (22) 
cos 2 sin 4 (xx) . 

Note that the wave vector of the A lg mode is completely fixed by and |k| (or E L ) as q x = — 2|k| cos 0. 
Phonons with different momenta are distinguishable in principle. We have different final states for 
different values, and to calculate the D band intensity, we need to sum over all possible final states by 
operating \Mji(Q)\ 2 with J Q 27r dO. Because J Q 27r cos 2 sin 4 OdO/ J Q 27r cos 6 OdO = 1/5, the polarization 
dependence of the D band intensity is written as 



where 9 in (9 out ) denotes the angle of an incident (scattered) electric field with respect to the armchair 
edge. 2 When the VV configuration (6* out = 9 in ) is used, the polarization dependence is approximated 
by I^ v (9 in ) ~ / cos 4 ^ n . On the other hand, when the VH configuration (^out = 9- m + n/2) is used, the 
polarization dependence is approximated by I^ H (9- m ) ~ / cos 2 6* in sin 2 9 in . These results are consistent 
with the experiment reported by Cancado et al [7]. Without a polarizer for the scattered light, we have 



2tt \ / r-2-K 

4qL„„2/Z1 \ , / / JQ„„„2n • 2nl„;„ 2/ 



a 



I D (9 in ) oc / rf0cos 4 cos 2 (# in ) + / dQ cos 2 © sin 2 © sin 2 (# in ). (24) 



Because J Q 2?r cos 2 © sin 2 OdO/ J Q 27r cos 4 OdO = 1/3, it may be rewritten as Ioi^m) oc cos 2 (# in ) + 

,2/ 



(1/3) sin (#i n ). With (Without) a polarizer for scattered light, the depolarization ratio Jd(90°)/ Id(0°) is 
1/5 (1/3). Generally, the D band polarization dependence is fitted using an empirical formula 

I D {9 iQ ) oc cos 2 (# in ) + bsm 2 {9 in ) + c, (25) 

where c is a constant, which probably originates from a defect beside the edge. The constants b and 
c determine the depolarization ratio: I D (90°) / Id(0°) = (b + c)/(l + c). When c is much larger than 
unity, the polarization behavior is obscured. If c is negligible, a smaller depolarization ratio (b = 1/5) is 
expected when a polarizer is used for the scattered light. 

4. Predictions of our model 

We have seen that the D band has a direct relationship to the bonding and antibonding orbitals (that is, 
an Ai g mode is excited through the first-order Raman process only from these orbitals). This conclusion 
has been derived based on two factors: the dominance of intervalley backward scattering (Sec. 3.1), and 



2 In calculating the D band Raman intensity, it is incorrect to sum over intermediate states specified by the electron's 
wave vector k, such as | ^ k Mji(9(k))| 2 . Because the phonon wave vector relates to the electron wave vector through 
q = —2k, actually means a summation for different phonons (final states). Thus, the Raman intensity is proportional to 

EJ^(e(q))| 2 or/ 2 ^e|M Jl (e)| 2 . 
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BZH by the armchair edge (Sec. 3.2). In this section, we report some consequences that are derived from 
these factors. 



4.1. The origin of Dispersive Behavior 



The D band frequency ud increases linearly with increasing excitation energy El as duo/dEL ~ 50 
cm~VeV, which is known as dispersive behavior [9,25-27]. In a previous paper [28], we pointed out 
that the dispersive behavior is mainly attributed to a quantum mechanical correction (self-energy) to the 
Ai 9 frequency. The modified energy of the A± g mode is written as ftw + ReU^(q,uj), where to is the 
bare frequency and 11^ (q, to) represents the self-energies of the A lg modes induced by electron-phonon 
interaction. The self-energy of an A\ g mode with q — |q| is defined as 



hw + shv^k — s'hv-p\\s. + q| + ie' 



(26) 



where e is a positive infinitesimal, /£ = lim /3 ^ OC) (l + e^^W ^) 1 is the Fermi distribution function 
with finite doping /i, and g s = 2 represents spin degeneracy. The electron-phonon matrix element 
squared \M\ 2 is constructed in Eq. (10) as \M\ 2 = ?f {1 - ss' cos(6'(k + q) - 6(k))}. 

In the continuum limit of k, II M (g, oj) is calculated analytically. The expression of the real part is 
given by 



ReYl^(q,uj)/G ep = Anfi 



hny u 2 — Vpq 2 9 t 



u)—v-pq 



+ hn\/v F q 2 



hvpq 



Hco + 2/i 

g ( ) + O hw-hv ¥ q 



Hco — 2a 
9 \ — I + 



sin 



1 i ftu) + 2/i 



hoj-\-hvpq Cf 



hv-pq 



\~ hvp-q-\-hoj 



sin 



2/i — hw 
hv-pq 



hv-pq 
(27) 
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Figure 6. A 3d plot of ReIl M (g, ui D ). The variables hvpq and \i are given in eV. Note that 
ReIl A1 (g, to) does not include the q dependence of the bare frequency. 




Hardening 

/ 47T/X 



1 .0 



1.0 



0.0 



where 9 X denotes a step function satisfying 9 X > = 1 and 6 X<0 = 0, g(x) = m(x + \fx 2 — 1), and G cp = 
g 2 p V/(2%hvF) 2 is a (dimensionless) coupling constant. A 3d plot of ReU fl (q,UD)/G ep is shown in 
Fig. 6. Interestingly, ReIl M (g, w D ) increases as we increase q. In fact, near the Dirac point, Rell /i (g, w D ) 
follows 



Since hvpq ~ the self-energy contributes to the dispersive behavior of the D (or 2D) band [9,25- 
27,29]. If we use G cp = 5 cm -1 , which is obtained from the broadening data published by Chen et 
al. [3], the self-energy can account for ~60% of the dispersion because ReU fl ^, (q,u q ) ~ G ep % 2 E L 
and G cp tt 2 ~ 30 cm -1 . It should be emphasized that the self-energy is calculated using only Eq. (10), 
and any artificial assumption, such as an adiabatic approximation, is not employed when calculating 
the self-energy. The physical origin of the dispersive behavior is easy to be understood with shifted 
Dirac cones [28]. In perturbation theory, the mechanism of dispersive behavior is almost the same as the 
mechanism where the G band exhibits hardening with increasing \fj,\. 

The dispersive behavior is not an inherent property of the D band but rather is a property of the Ai g 
mode. The armchair edge is involved in the activation of the D band. However, the dispersive behavior 
itself has nothing to do with the edge. In other words, there are processes that can exhibit dispersive 
behavior, besides the Raman D band. A good example is the 2D band, for which dispersive behavior is 
observed because the 2D band consists of two A lg modes. 

4.2. Intravalley Phonons 

In addition to the zone-boundary A\ g mode, we can determine the effect of BZH on zone-center 
(intravalley) optical phonon modes: BZH forbids an intravalley transverse optical (TO) mode to appear 



Reil^o(<?,w) ~ G cp 7r 2 hv F q. 



(28) 
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as a prominent Raman band at the armchair edge. The Dirac equation is the most useful way of showing 
this. The electron-phonon interactions for the LO and TO modes with (nonzero) momentum (q X) q y ) are 
written as 



3 iqr 



>»—[ n " ll:i (] rwi " _ : _ ) (29) 




# TO (q) = Scp^ |^ + Q ° V% __"__), (30) 

and their matrix elements are obtained from Eq. (17) as 

M LO (q) = y^5(-k' y + q y + k y )x 

{S(-k' x + q x + k x ) \{a x )q y - (<7 y )q x \ + S(-k x - q x + k x ) \{cr x )q y + (<J y )q x \} , (31) 



M TO ( q ) = ^5(-k' y + q y + k 



IX 



|q| 

+ Qx + k x ) [(a x )q x + (<7 y )g v ] + S(-k x - q x + k x ) [(a x )q x - (a y )q y }} , (32) 

where (a x ) = (se ie> + s'e~ e )/2 and (a y ) = -i(se ie ' - s'e-°)/2. For the TO mode, M TO (q) vanishes 
in the q x — > limit, due to the interference between the valleys. In the q y — > limit, 6' = ix — holds, 
and the orbital of the matrix element, (a x ), vanishes when s' = s. Thus, the electron-phonon matrix 
element for the T point TO mode is suppressed compared with that of the LO mode: the T point TO 
mode is missing in the G band at the armchair edge [10,20,30,31]. 

4.3. D Band Splitting 

The D band is composed of the (two) A\ g modes that are emitted from two electronic states with 
bonding or antibonding orbitals (0 = or it). This fact results in the splitting of the D band if the 
trigonal warping effect is taken into account [12]. The splitting width increases with increasing incident 
laser energy E L as 

Au D = 25 [cm" 1 ], (33) 

where 7 (~ 3 eV) is the hopping integral between nearest neighbor atoms. This formula is derived by 
noting that the wave vectors for the two A lg modes that originate from the bonding and antibonding 
orbitals, are different due to the trigonal warping effect. The difference is estimated with the lattice 
tight-binding model as [12] 

A^(^V. (34) 

a V37/ 

For simplicity, let us assume that the phonon dispersion relation for the D band is isotropic about the 
Dirac point. From Fig. 7 it is clear that the two phonon modes with q and q n have different phonon 
energies, which results in the double peak structure of the D band. The difference between the energy 
of the phonon mode with q and that with q n is approximated by 

Au D = ^Aq, (35) 
oq 
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where dup/dq is the slope of the phonon energy dispersion. We interpret the dispersive behavior that 
occurs as a result of the dispersion relation of the phonon mode. Then we have 

du D _ dE L du D 
dq dq dE L 

Putting \dEi,/dq\ = Hvp in Eq. (36), and combining it with Eqs. (34) and Eq. (35), we obtain Eq. (33). 

Figure 7. Two principal wave vectors of the phonons (q and q n ) contributing to the D band. 
The phonon dispersion relation near the Dirac point. The energy difference between cop(q ) 
and u D (q 7T ) appears as two peaks in the D band. 




5. Prospects 

The D band has a close relationship with the energy gap. A direct relationship between the zone- 
boundary Ai g mode and the energy gap is seen in the Pierls instability at armchair nanoribbons [32,33]. 
In the Dirac equation, an energy gap is represented by a Dirac mass term that appears in the off-diagonal 
terms as 

fj _ ( ' P m<T * 
\ m*a x Vpcr' ■ p 

The spectrum has an energy gap 2\m\, because the energy dispersion is given by ±*/ (hvpk) 2 + \m\ 2 . 
The zone-boundary A lg mode is related to the Dirac mass because the electron-phonon interaction 
[Eq. (11)] appears as a mass term in the Dirac equation: m — > g op e~ 4q r . The dominance of intervalley 
backward scattering becomes clear with respect to the mass term, because a nonzero mass tends to stop 
a massless particle by backward scattering. 

Interestingly, the armchair edge itself can be modeled as a singular mass term in the Dirac equation, 
whereby Eq. (17) is obtained as a solution of the model [34]. This theoretical framework leads us to find 
a relationship between the dominance of intervalley backward scattering represented by Eq. (10) and 
BZH. To see a connection between them, it is important to recognize that in deriving Eq. (10), the orbital 
structures of Eqs. (1) and (2) play a decisive role. These orbitals have to relate to each other by mirror 
symmetry with respect to the x-axis [x — > —x or k x — > —k x ] in Fig. 1. In fact, Eq. (2) is constructed by 
replacing 6 with n — 0' in Eq. (1). The orbital should be invariant under this replacement. 
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Defects, such as a lattice vacancy and a topological defect, are considered to be sources that 
increase the D band intensity [c-term in Eq. (25)]. This speculation is reasonable, because creating 
a lattice vacancy inevitably involves the antibonding orbital. Unfortunately, the wave functions and 
the corresponding BZ in the presence of defects are difficult to construct rigidly in the framework of 
a tight-binding lattice model. This difficulty prevents us from calculating the electron-phonon matrix 
element exactly. However, there is the possibility of obtaining a good approximation of the matrix 
element using the Dirac equation. 

6. Conclusion 

An orbital and wave vector are the basic idea for a molecular and a crystal, respectively. Knowing 
the correlation between them is the key to understand the D band. The D band originates from two 
orbitals: the bonding and antibonding orbitals on the A^-axis. Observing the D band is the same thing as 
selecting the two orbitals from the various orbitals that compose the iso-energy section of the Dirac cone. 
This idea leads us to expect the optical control of edge chiralities [35]. A slight asymmetry between the 
bonding and antibonding orbitals, induced by the trigonal warping effect, may be important in terms of 
understanding the stability of the armchair edge under laser light irradiation. A closer study of the D 
band based on the Dirac equation will be fruitful for a further investigation of the physics of the D band. 

Acknowledgments 
References 

1. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Katsnelson, M.I.; Grigorieva, I.V.; 
Dubonos, S.V.; Firsov, A. A. Two-dimensional gas of massless Dirac fermions in graphene. 
Nature 2005, 438, 197. 

2. Zhang, Y.; Tan, Y.W.; Stormer, H.; Kim, P. Experimental observation of the quantum Hall effect 
and Berry's phase in graphene. Nature 2005, 438, 201. 

3. Chen, C.F.; Park, C.H.; Boudouris, B.W.; Horng, J.; Geng, B.; Girit, C; Zettl, A.; Crommie, 
M.F.; Segalman, R.A.; Louie, S.G.; Wang, F. Controlling inelastic light scattering quantum 
pathways in graphene. Nature 2011, 471, 611 . 

4. Ando, T; Nakanishi, T; Saito, R. Berry's Phase and Absence of Back Scattering in Carbon 
Nanotubes. J. Phys. Soc. Jpn. 1998, 67, 2857. 

5. Tuinstra, F.; Koenig, J.L. Raman Spectrum of Graphite. J. Chem. Phys. 1970, 53, 1126-1130. 

6. Katagiri, G.; Ishida, H.; Ishitani, A. Raman spectra of graphite edge planes. Carbon 1988, 
26, 565-571. 

7. Cancado, L.G.; Pimenta, M.A.; Neves, B.R.A.; Dantas, M.S.S.; Jorio, A. Influence of the Atomic 
Structure on the Raman Spectra of Graphite Edges. Phys. Rev. Lett. 2004, 93, 247401. 

8. You, Y; Ni, Z.; Yu, T; Shen, Z. Edge chirality determination of graphene by Raman 
spectroscopy. Appl. Phys. Lett. 2008, 93, 163112. 

9. Gupta, A.K.; Russin, T.J.; Gutierrez, H.R.; Eklund, PC. Probing Graphene Edges via Raman 
Scattering. ACS Nano 2009, 3, 45-52. 



Version November 22, 2012 submitted to Crystals 



17 of 18 



10. Cong, C; Yu, T.; Wang, H. Raman Study on the G Mode of Graphene for Determination of Edge 
Orientation. ACS Nano 2010, 4, 3175-3180. 

11. Thomsen, C.; Reich, S. Double Resonant Raman Scattering in Graphite. Phys. Rev. Lett. 2000, 

85, 5214-5217. 

12. Sasaki, K.; Kato, K.; Tokura, Y.; Suzuki, S.; Sogawa, T. Pseudospin for Raman D band in 
armchair graphene nanoribbons. Phys. Rev. B 2012, 85, 075437. 

13. Wallace, P.R. The Band Theory of Graphite. Phys. Rev. 1947, 71, 622-634. 

14. Slonczewski, J.C.; Weiss, P.R. Band Structure of Graphite. Phys. Rev. 1958, 109, 272-279. 

15. Sakurai, J. Advanced Quantum Mechanics; Addison- Wesley: Canada, 1967. 

16. Baranov, A.V.; Bekhterev, A.N.; Bobovich, Y.S.; Petrov, V.I. Interpretation of certain 
characteristics in Raman spectra of graphite and glassy carbon. Optics and Spectroscopy 1987, 
62, 612-616. 

17. Pocsik, I.; Hundhausen, M.; Koos, M.; Ley, L. Origin of the D peak in the Raman spectrum of 
microcrystalline graphite. Journal of Non-Crystalline Solids 1998, 227230, Part 2, 1083 - 1086. 

18. Sasaki, K.; Saito, R. Pseudospin and Deformation-Induced Gauge Field in Graphene. Prog. 
Theor. Phys. Suppl. 2008, 176, 253-278. 

19. Sasaki, K.; Wakabayashi, K.; Enoki, T. Electron Wave Function in Armchair Graphene 
Nanoribbons. J. Phys. Soc. Jpn. 2011, 80, 044710. 

20. Sasaki, K.; Saito, R.; Wakabayashi, K.; Enoki, T. Identifying the Orientation of Edge of Graphene 
using G Band Raman Spectra. J. Phys. Soc. Jpn. 2010, 79, 044603. 

21. Maeta, H.; Sato, Y. Raman spectra of neutron-irradiated pyrolytic graphite. Solid State Commun. 
1977, 23, 23 - 25. 

22. Nemanich, R.J.; Solin, S.A. First- and second-order Raman scattering from finite-size crystals of 
graphite. Phys. Rev. B 1979, 20, 392-401. 

23. Griineis, A.; Saito, R.; Samsonidze, G.G.; Kimura, T; Pimenta, M.A.; Jorio, A.; Filho, A.G.S.; 
Dresselhaus, G.; Dresselhaus, M.S. Inhomogeneous optical absorption around the K point in 
graphite and carbon nanotubes. Phys. Rev. B 2003, 67, 165402. 

24. Sasaki, K.; Kato, K.; Tokura, Y; Oguri, K.; Sogawa, T. Theory of optical transitions in graphene 
nanoribbons. Phys. Rev. B 2011, 84, 085458. 

25. Vidano, R.; Fischbach, D.; Willis, L.; Loehr, T. Observation of Raman band shifting with 
excitation wavelength for carbons and graphites. Solid State Commun. 1981, 39, 341 - 344. 

26. Matthews, M.J.; Pimenta, M.A.; Dresselhaus, G.; Dresselhaus, M.S.; Endo, M. Origin 
of dispersive effects of the Raman D band in carbon materials. Phys. Rev. B 1999, 
59, R6585-R6588. 

27. Casiraghi, C; Hartschuh, A.; Qian, H.; Piscanec, S.; Georgi, C; Fasoli, A.; Novoselov, K.S.; 
Basko, D.M.; Ferrari, A.C. Raman Spectroscopy of Graphene Edges. Nano Lett. 2009, 9, 1433. 

28. Sasaki, K.L; Kato, K; Tokura, Y; Suzuki, S.; Sogawa, T. Decay and frequency shift of both 
intervalley and intravalley phonons in graphene: Dirac-cone migration. Phys. Rev. B 2012, 

86, 201403. 

29. Piscanec, S.; Lazzeri, M.; Mauri, F; Ferrari, A.C; Robertson, J. Kohn Anomalies and Electron- 
Phonon Interactions in Graphite. Phys. Rev. Lett. 2004, 93, 185503. 



Version November 22, 2012 submitted to Crystals 



18 of 18 



30. Begliarbekov, M.; Sul, O.; Kalliakos, S.; Yang, E.H.; Strauf, S. Determination of edge purity in 
bilayer graphene using mu-Raman spectroscopy. Appl. Phys. Lett. 2010, 97, 031908. 

31. Zhang, W.; Li, L.J. Observation of Phonon Anomaly at the Armchair Edge of Single-Layer 
Graphene in Air. ACS Nano 2011, 5, 3347-3353. 

32. Fujita, M.; Igami, M.; Nakada, K. Lattice Distortion in Nanographite Ribbons. Journal of the 
Physical Society of Japan 1997, 66, 1864-1867. 

33. Son, Y.W.; Cohen, M.L.; Louie, S.G. Energy Gaps in Graphene Nanoribbons. Phys. Rev. Lett. 
2006, 97, 216803. 

34. Sasaki, K.; Wakabayashi, K. Chiral gauge theory for the graphene edge. Phys. Rev. B 2010, 
82, 035421. 

35. Begliarbekov, M.; Sasaki, K.I.; Sul, O.; Yang, E.H.; Strauf, S. Optical Control of Edge Chirality 
in Graphene. Nano Letters 2011, 11, 4874-4878. 

© November 22, 2012 by the authors; submitted to Crystals for possible open access 
publication under the terms and conditions of the Creative Commons Attribution license 
http://creativecommons.Org/licenses/by/3.0/. 



